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ABSTRACT 

Parameter  plane  techniques  are  applied  to  compensation  of  third 
order  control  systems  by  extending  the  principle  to  include  the  quad- 
ratic form  of  the  parameters. 

Equations  for  constant  bandwidth  curves  on  the  parameter  plane 
are  developed. 

Recursive  equations  for  the  extension  of  these  methods  to  higher 
ordered  systems  are  presented. 

A  computer  program  for  the  computation  and  plotting  of  constant 


zeta,  omega  and  bandwidth  loci  is  written. 
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1 .  Introduction 

Until  1959  when  Mitrovic  [lj  introduced  what  is  now  called  the 
Coefficient  Plane,  analysis  and  compensation  of  control  systems  were 
accomplished  primarily  by  frequency  response  and  root  locus  tech- 
niques. With  the  advent  of  the  Coefficient  Plane  there  was  then  a 
method  available  which  allowed  the  control  engineer  to  determine  the 

root  locations  of  a  plant  as  two  coefficients  of  the  characteristic 
I 

equation  were  varied. 

At  about  the  same  time  E.K.  Ross  and  T.C.  Warren  [_2j  in 
conjunction  with  Dr.  G.J.  Thaler  published  a  numerical  method  of 
determining  single  and  multiple  section  cascade  compensation  while 
maintaining  a  constant  error  coefficient.  Shortly  thereafter  in 
i960  CD.  Pollak  [3J  and  Dr.  Thaler  published  yet  a  different 
numerical  method  for  stating  compensator  pole  locations  for  a 
given  gain. 

In  1964>  Siljak  [4  J  extended  Mitrovic fs  work  to  apply  to  pro- 
blems in  phich  the  system  variables  or  parameters  appeared  as  a 
product  in  one  or  more  coefficients.  This  paper  will  further  extend 
the  Parameter  Plane  concept  to  include  the  quadratic  form  of  the 
coefficients.   Specifically  the  case  of: 

S  a^k  =  0  (1-1) 

k-0 


where: 


and: 


•k  ■  \°<*  +  v* +  vy3  +  v3  +  y3* +  %      (i-2) 


b,  ,c,  ,<L  ,e,  ,f,  ,g,   are  real  and  constant 
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o<  and  A  are  real  and  variable 

These  curves  which  will  be  plotted  on  the  ex   -  A      or  Parameter 
Plane  will  specify  two  of  the  roots  of  Equation  (l-l)  for  any 
choice  of   CX   and  A   . 
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2.  Concepts  of  the  Parameter  Plane 

The  Parameter  Plane  was  derived  by  a  study  of  the  character- 
istic equations  of  control  systems  where  two  parameters  were  allowed 
to  vary.  Here  a  presentation  of  the  basic  ideas  of  Siljak  ^4j  will 
be  studied. 

Consider  a  characteristic  equation  of  the  form: 

i:  vk  -  °  (") 

k=0  K 

where  the  coefficients  a,,  a0,...,  a  are  real  and  s  is  the 

l7  d7        7  n 

Laplace  Transform  variable: 

s  =  CT  +  jo)  -  -  3^0  +  o  uS]   1-  J2  (2-2) 

By  forming  powers  of  s  it  is  shown  \_A \  that: 

Sk  -  60k  TK(-3)  *  *V  1-  J2  '  Uk(-3)  (2-3) 

where : 

V-5)  =   (-i)\(5) 

(2-4) 

\(-S)  =   (-Dk+1uk(3) 

The  U,  and  T,  are  the  first  and  second  kind  of  Chebyshev 
functions  and  are  generated  by  the  relations  [4I  : 


T 


k+1 


(3)  -25Tk(3)  +  Vi<3>    ■  ° 


with: 


uk+1(3)  -23uk(3)  ♦  ^(Z)    =  0 

t0(3)    -   1  uo(3)    =   0 

tx(3)    =  0  Pj(3)   -  -1 


(2-5) 
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Taking  Equation  (2-1),  substituting  Equation  (2-3)  in  it,  then 
equating  the  real  and  imaginary  parts  to  zero  separately  and  employ- 
ing Equations  (2-4)  and  (2-5)  we  will  obtain  the  expressions  in 
terms  of  Chebyshev  functions: 

,k    k.. 


i:  {-l)\u>\A3)      -  0  (2-6) 

k=0       lc    k  X 

E   (-l)\co\(3)   =  0  (2-7) 

k=0       *    K 

Now,  assuming  as  Siljak  did,  that  the  a,   coefficients  of 

Equation  (2-l)  are  linear  functions  of  o<      and  /3  : 

■t  ■>*  +0k/+ak  <2-8> 

with: 

b,  ,c,  ,d,    constant  and  real 

cK   and  y^    variable  and  real 
we  may  then  substitute  Equation  (2-8)  into  Equations  (2-6)  and 
(2-7)  and  obtain: 

L/3    =-: 

+  C2/3     = 

where : 


V  +cifi  °-\ 


h~  +  V3  =  _B2 


(2-9) 


b,  =  s  c-i)\« v ,  b ....  t:  (_i)\«*>v 

1       k=0               K         k  -1  ^       k=0               k        E 

B,   -  £    (-l)\^kUk.1  B2   .  p    (-D  V\ 

k=0  k=0 
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Now  Equations  (2-9)    may  be   solved  simultaneously  for    crt     and  A 
giving: 

ClD2  -  CA  B2D1  -  B1B2         ,211, 

"  B1C2  "  32C1       I        =  D1C2  "  B2C1 

Equations  (2-11)  now  give  crt  and  /^  in  terms  of   o   and  60  and 
either  may  be  fixed  and  the  other  varied  from  zero  to  infinity  which 
will  give  a  curve  on  the  o<  -  /3  plane.  This  may  be  repeated  for 
as  many  values  as  needed  giving  an  entire  family  of  curves  on  the 
ex  — /S     plane  which  completely  describe  the  behavior  of  a  pair 
of  complex  roots.  These  curves  will  usually  be  constant  zeta  or 
constant  omega  loci. 

If  s  =  C3"  is  substituted  in  Equation  (2-1)  and  the  same 
procedure  is  followed,  curves  of  constant  sigma  or  values  of  s 
along  the  real  axis  may  be  obtained. 

This  then  is  the  general  procedure  to  be  followed  in  the  devel- 
opment which  follows  for  the  quadratic  case. 
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3.  Cascade  Compensation  of  Third  Order  Systems 

The  initial  investigation  of  the  ex  - p    quadratic  case  is  made 
by  attempting  to  compensate  a  specific  third  order  system  with  two 
identical  sections  of  cascade  compensators  maintaining  a  constant 
velocity  coefficient.   The  uncompensated  plant,  Figure  3-1  >  is  chosen 
so  as  to  afford  not  only  lead,  hut  lag  compensation,  and  with  this 
selection,  the  system  gain  is  set  to  place  the  complex  roots  at  the 
stability  limit,  which  also  insures  the  need  for  double  section  comp- 
ensation over  a  large  range  of  zeta  and  omega.  Establishing  the 
generalized  compensated  plant,  Figure  3-2,  the  characteristic 

equation  is  determined: 

2  5        tc  2       o     2\   4        /c   2        no     2n    3        /t^     2 
z   s^  +   (6z     +  2pz   )s     +   (5z     +  12pz   Js     +   (lOpz     + 

6p2z2  +  20p2)s2  +   (5p2z2  +   60p2z)s  +  30p2z2     =0  (3-1) 

Simplifying  this  equation  and  making  a  substitution  of  variables: 

/3    =     p  <X     =     p/z  (3-2) 

an  equation  in  terms  of    <X    and     yo      is  obtained: 

s5   +   (6  +  2/3)s4  +   (5     +  12/3   +/32)s3  +   (10/3  +  6/32  + 

30o<2)s2  +  (5/3  2  +  30oyg)s  +  30/52     =     0  (3-3) 

and  when  the  terms  are  arranged  in  terms  of  the  parameters: 

c*2(30s2)   +    oy3(60s)   +  y3(2s4  +  12s3  +  10s2)    + 

/32(s3  +  6s2  +  5s  +  30)   +  (s5  +  6s4  +  5s3)      =     0  (3-4) 

which  shows  the  general  coefficient  form: 

he*/5  +  ek/3  -  fy3' +  sk) 


S    (bk**  +  cka  +  &v<*/2>    +  ev/3  +  fv/6     +  g*)      =0       (3-5) 
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UNCOMPENSATED     SYSTEM 
FIG.    3-! 


-±o 


-4 


Gc 


P     (S^2) 


30 


S(S*I)(S*5) 


COMPENSATED       SYSTE  M 
FIG^  3-2 
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or: 

n 


k 


U   VT    =    o  (3-6) 

k=0 

■where  all  of  the  cn  's  are  zero, 

k 

Employing  Equations  (2-6)  and  (2-7) 

i:  (-D\^\^(t)  -  o  (2-6) 

k=0       k         k  l 

jt   (-DV\(2)   =  0  (2-7) 

k=0       *    * 

Equation  (3-4)  may  be  split  into  a  real  and  imaginary  part  in  terms 

of  3  and  60  : 

^2b21  +  0^/3  \  +/3e11  +/32f21  +  g1    =    0  (3-7)1 

o<2B22  +    cx/5D2   +  /3E12  +  /32F22  +  G2     =     0  (3-8) 


where : 

B21   =      30^^(3) 

B22  =     30oi2U2(J) 

D1     --60coUo(3) 

L2     =  -60*0^(1) 

(3-9) 

(3-10) 

(3-11) 
(3-12) 
B11     =  2  604U3('S)   -  12dJ3U2(3)    +  10  0)^(3)  (3-13) 

E12  =  2cJlJ4C$)   -  12^o3U3("S)    +  10o32U2(f  )  (3-14) 


The   second  number  on  the  double   subscripted  variables  and 
the  number  of  the   single   subscripted  variables  indicates  the  real(l) 
or  imaginary (2)   part  of  the  equation. 
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F21   -  -«o3U2(S)    +  6o)2U1(  S)   -  5o3Uo(S)   +30U_1(J)  (3-15) 

F22  "  -^3u3(^)   +  6c°2u2(  ^  ~  5^UX<J)   +  30UQ  (J)  (3-16) 

Gl     =-^5u4(3)  +  6aAj3(3)  -56o\-2(S)  (3-17) 

G2     =  -oi5U5(S)   +  6«A4(J)  -  5^3U3(J)  (3-18) 

Further  collecting  the  terms  of  Equations  (3-7)  and  (3-8)  we 
obtain: 

<*2b21  +c*Pi  +  Qi     =     °  (3-19) 

<*  2BOQ  +  <X  P0  +  Q0     =     0  (3-20) 


where : 


22         "    2        ^2 


P-L     -/^D]L  (&-21).; 


P2     =/3D2  (3-22) 

S.  'Z3  En  +/^2p2i +  Gl  (3~23) 

Q2  =/3E12  +/32f22  +  G2  (3~24) 

We  now  have  obtained  a  set  of  two  non-linear  algebraic  equations  in 
terms  of    o     and   ex)    which  must  be   solved  simultaneously  for  the 
correct  values  of    t*    and  A    .     To  do  this,    the  method  with  the  best 
chance  of  sucess  appears  to  be  Sylvester's  Method  in  which  we  form 
a  set  of  four  equations  by  taking  the  original  Equations  (3-19)   and 
(3-20)   and  forming  two  more  by  a  multiplication  with  o{    giving: 

CX2B21  +cXP1  +  0^     =     0  (3-19) 

CX2B22  +<*P2  +  Qg     =     0  (3-20) 
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°<  B21  +  <*  Pl  +  ^Ql   =  ° 


o<3B22  +  CX2F2   +   o(Q2   =  0 


Now  placing  these  equations  in  matrix  form: 


0    B     P    0, 


0    B£2   P2   Q2 


B21  Pl 


Qj   0 


B22  P2    Q2   0 


r^3 


=  0 


we  may  say  that  if  the  ck  '  s  are  not  zero  then: 


B 


21 


Qn 


0 

B22 

P2 

Q 

B21 

P! 

«1 

0 

B22 

P2 

Q2 

0 

0 


(3-25) 
(3-26) 


(3-27) 


(3-28) 


Expanding  this  determinant  "we  obtain: 

~B21Q2  +  B21B22Q1Q2  +  P1P2B21Q2  "  P1Q2B22  +  Q1Q2B21B22 
-Q1B22  -  Q1P2B21  +  SP1P2B22  =  ° 


(3-29) 


Substituting  Equations  (3-21)  through  (3-24)  in  Equation  (3-29)  we 
achieve  a  fourth  order  equation  in  /3   : 
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pH-4 


2  ,2 

321  +  2P21P22B21B22  +  V2P22B21  -  P22B1B22  " 

P21B22  "  -?21I,2B21  +  D1B2P21B22)   + 


^3(-2B12F22Bll  +  2E11P22B21B22  +  2P21E12B21B22  + 


B1B2E12B21  "  B1E12B22  "  2E11P21B22  " 


D2E11B21  +  WllB22>  + 


fi   (-E12B21  -  2F22G2?21  +  2EnE12B21B22  J  ZP^B^S^   + 
2G1P22B21B22  +  B1B2G2B21  "  D1G2B22  "  E11B22  " 
2P21°1B22  -  D2G1B21  +  WlB22)   + 
/3  (-^B^Va  +  2E11°2B21B22  +  2G1E12B21B22  ~  2E11G1B22>   + 

("G2B21  +  2G1G2B21B22  -  °M>>  *  °  O30) 

from  -which  the  coefficients  may  be  determined  by  a  substitution  of 
;^ttalffes  (3-9)  through  (3-18)  and  the  values  of  the  first  kind  of 
Chebyshev  functions,  Equation  (2-5)  i  in  terms  of  o  and  cO  .  Since 
the  solution  of  a  fourth  order  equation  is  at  best  difficult,  it  is 
at  this  point  a  digital  computer  becomes  a  necessity. 

The  major  problem  is  not  the  actual  solution  of  the  quartic 
itself,  but  rather  the  proper  choice  of  one  of  the  four  solutions. 
There  are  two  marked  characteristics,  however,  which  help  in  the 
selection.  These  ares 

a.  Complex  answers  to  the  quartic  have  no  physical  significance 
and  may  therefore  be  discarded  as  erroneous. 
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b.  The  definition  of  c*   requires  that  <*  and  /3    be  of 

the  same  sign  so  that  p  and  z  will  be  of  identical  sign. 
Using  this  information  and  that  available  from  the  Ross-Warren  £2  J 
method  (Appendix  I)  as  to  compensator  pole  and  zero  location,  it  is 
found  that  the  solution  to  the  /3    quartic  (Appendix  II)  is  the 
largest,  positive,  real  value. 

liaw   entering  Euqation  (3-21)  with  this  value,  and  evaluating  the 
other  coefficients  (Appendix  II )  we  find  that: 


<*  = 


[-  Qi  /  B21]  i  (3-31) 


for  in  the  third  order  case  P.,  is  always  identically  zero. 

Thus,  with  the  programing  of  the  appropriate  equations,  the 
digital  computer  could  give  all  of  the  values  and  plot  the  constant 
zeta  and  constant  omega  loci  on  the  Parameter  Plane  for  any  desired 
values.  With  the  ex.   and  /o    values,  Figure  3-3,  and  the  graphs  of 
the  Parameter  Plane  Lag,  Figure  3-4>  and  Lead,  Figure  3-5,  compen- 
sation zones,  a  spot  check  of  the  velocity  coefficient  is  made  with: 

where  K  is  the  compensated  root  locus  gain. 
c 

Case  1.     3  =  .5  ua  =  1.0     r  =  -.5  +  j  0.86603 

p  =  /3  =  0.03037  c*  =  p  /zn  =  0.36895   z  =  0.08231 


(1.0) (.52  +  .866032)^(4>52  ±   .866032)*( .469632  +  .866032) 


(.417692  +  .866032) 


K  =  4.81100 
c 


K  =  (4.81100) ( .08231 r    m      7>o6 
(1.0)(5.0)(.03037)2 
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THE  ALFA  VALUES  ARE 


ZETA      .000000     .100000     .200000     .300000     .400000 
..50_.  .23818E*00  _.25986E*00  . 27452E*00. . 28300E*00  ,28584E*00 


1.00  .43121E*00  .45283£*00  .45547E*00  .44130E*00  .411266*00 

1.50  .64650E*00  .65840E*00  .62918E-00  .58083E*00  .51866E*00 

2.0  0  .8819.7E  +  0  0.  .  857 Q6E*0  0_.. 80  0  0 4E*0  0__.  72553E*0  0.  63375E*0JL 

2.50  .11372E*01  .13493E*01  .22198E+01  .47367E*01  ,1Q206E*Q2 


3.00  .14126E*Q1  ,15121E*Q1  .17341E-01  ,20955E*Q1  .25816E*01 
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Case  2.   S  =  .2    cO   =  3.0    r  =  -.6  +  j  2.93943 

/3  =  Pc  =  4.5954     *  .  pc  /  zc  =  1.7341  zc  -  2.65002 

K  (3.0H.42  t   2.939432)*(4.42  +  2.939432)'ir(3.99542  ±   2.93943fe) 
C  (2.050022  +  2.939432) 

K  =  90.21638 
c 

K  ■  (^>21638)(2.65002)2   =  6#00021 

v   (l.O)(5.0)(4.5954) 

Case  3>    3  -  .8     <*   =  4.5    *  -  -3.6  +  j  2.7 

/3  =  Pc  =  10.191   <*  =  Pc  /  zc  =  3.373   zc  =  3.02134 

K  (4.5) (2. 62  +  2.72)*(l.42  +  2.72)*(6.5912  *   2.72) 
c         (.578662  +  2.72) 


K,  =  341.3239 


2 


*■  (341.323?) (3.02134)^   H  6#00015 
(l.o)(5.o)(io.l9J) 


Noting  that  all  of  the  calculated  K * s  are  greater  than  or  equal  to 
6,0,  the  velocity  coefficient  of  the  uncompensated  system,  a  check 
of  the  root  loci  is  made.  As  shown  in  Figures  3-6,  3-7  and  3-8, 
the  desired  root  location  is  also  achieved. 
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4.  Generalization  of  the  Third  Order  System 

In  Section  3,  the  development  of  a  specific  type  one,  third 
order  system  was  shown.  Here  a  generalized  development  of  the 
equations  necessary  to  write  a  computer  program  for  any  third 
order  system  is  shown. 

Starting  as  before  by  specifying  the  compensated  plant,  (this 
time  however  in  generalized  terms),  as  in  Figure  4-1,  the  character- 
istic equation  is  written. 

s5  +  (2p  -i-Er^s4  +  (p2  +  2VZt±  +S7Tri)s3  + 

(p2Zr.    +  SpETTr.    +1Tr.)s2  +  (p^TPr.    +  2pTTr.)s     + 
1  2.     1  l  £     l  l 

p^r.   +  K(2g-  s2  +  2-2-  8  +  p2)  (4_i) 

z 


where : 


Dr.  =  the  sun  of  the  roots 


i 


STPr.   =  the  sum  of  the  product  of  the  roots  taken 
"*  n  at  a  time 

TFr.   =  the  product  of  the  roots 

¥ith  the  substitution  of: 

/3   =  p         c*  =  p/z  (4-2) 

and  the  collection  of  like  terms  in  o<    and  A   ,  we  obtain: 

<*2(Ks2)    +  <x/3(2Ks)    +/32(a3  f^r.s2  +2TJVS  +TTx±  +  K) 

+  /3(2s4  +  2Zr.s3  +  2ETTr.s2  +  2TTr.s) 
/^v  1  2     1  1    ' 

+  (s5  +  Zris4+Zrn,ris3  +  TTris2)  (4-3) 
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Employing  Equations  (2-6)   and  (2-7)   and  defining: 
2T 


B21  =  X  o»\($)  (M) 

^o(S)  (4-5) 


B22   =KW.2 


D1     =  -2Kca)Uo(  $  )  (4-6) 

D2     --ZKcaU^S)  (4-7) 

En  =  2oAj3(S)   -22^03^(3)  +  2Z,niri^2U1(J) 

-  2lTri^U0(S)  (4-8) 

E12  =  2£°VS)  -eZ»i**S<3i  ♦  2LTr.602U2(j) 

-  2lT:ria>U1(,5)  (4-9) 

F2i  *-*\i3)  +Z  *!*>%(  5)  -rpi^u0(5) 

+  (TTr.    ♦  ^(J)  (4-10) 

F22  "  '^\(S)   +2^^2U2(S)  -ZTTr.coU^S) 

*  (Tr.    +K)Uo(S)  (4-11) 

<»1     =  -6o\(  3  )   *  2r.aA3(  !  )  -Il7Tric03U2(  $  ) 

-hTTr^U^J)  (4-12) 

G2     =  -^(1)   +S  riu3\(  1)  -S^r.^d  )  . 

+  ^^00^(1)  (4-lS) 

Px     =  /&\  (4-14) 

P2     =  /Sd2  (4-13) 

\     =   /3En  h-^F^  +  G1  (4-16) 

Q2     =    /^E12  +/^2f22  +  G2  (4"17) 
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we  find: 

<X2B21  +  o<P1  +  Q1  =  0  (4-18) 

C*2B22  +  <7VP2  +  Q2  =  0  (4-19) 

which  are  two  non-linear  algebraic  equations  completely  general- 
ized in  terms  of  the  uncompensated  system  poles  and^root  locus  gain, 

o  ,  £&     and  the  first  kind  of  Chebyshev  Functions.  Equations 
(3-30)  and  (3-31)  are  now  employed  to  obtain  the  solutions  for  o^ 
and  /o    ,  using  the  same  selection  rules  as  before. 

Several  third  order  systems  were  investigated  by  the  application 
of  the  generalized  equations  and  the  Parameter  Plane  curves,  Figures 
4-2  through  4-9  >  were  plotted.  Of  these,  the  K/s  family  appears 
the  most  interesting.   Further  investigation  of  three  of  the  curves 
in  this  family,  Figures  4-2,  4-3  and  4-4  shows  that  there  is  a  rela- 
tionship between  K,  the  root  locus  gain,  o<  and  3   • 
These  relations  are: 

a.  Choose  a  point  on  the  l/s   ex  —ft   plane. 

b.  Zeta  reads  directly. 

c.  Determine  the  actual  omega  at  that  point  by  multiplying 
the  value  read  by  the  cube  root  of  the  uncompensated 
system  gain. 

d.  Read  the  value  of  c*  directly  from  the  point  chosen. 

i 

e.  Eead  the  value  of  /3  from  the  point  chosen. 

f.  Obtain  the  true  value  of  p    by  multiplying  this  value  by 
the  cube  root  of  the  uncompensated  system  gain. 

By  this  method,  the  values  of  c\  and  /3  may  be  determined  for 
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all  K/s  systems  from  the  one  universal  curve. 

Although  the  majority  of  the  curves  are  easily  interpreted, 
there  are  some,  Figures  4-9  and  4-10,  which  require  careful  corre- 
lation of  the  data  points  and  the  curves.  As  shown  in  Figures  4-11 
and  4-12,  the  correct  curves  may  be  picked  out  and  redrawn  if 
necessary  for  further  clarity  and  ease  of  use. 
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5.  Bandwidth  Curves  on  the   o<  ->S  Plane 

In  many  instances,  there  is  also  a  bandwidth  criterion  imposed 
on  the  engineer  as  well  as  an  optimal  operating  point  for  the  plant 
under  consideration.  With  this  in  mind,  equations  for  the  plotting 
of  constant  bandwidth  curves  on  the  oL  - /S  plane  are  developed. 
For  the  purpose  of  this  development  a  constant  bandwidth  curve  will 
be  defined  as: 

A  constant  bandwidth  curve  for  Gf(j*o  )  =  M  is  a  curve  drawn 

upon  the  parameter  plane  which  specifies  the  relation  between 
the  parameters  necessary  if  the  transfer  function  G(s),  which 
is  a  function  of  the  parameters,  is  to  have  magnitude  M  at 
the  real  frequency  60  ^.  ^5~] 

Once  these  curves  are  obtained  they  may  be  superimposed  on  the 

parameter  plane  thus  indicating  what  values  of  the  parameters  are 

necessary  in  order  to  meet  the  specifications. 

Taking  the  rational  transfer  function  and  defining  it: 

m       m-1 
p/_\     P_s  -i-p  ,s    +  ...p  s  +  p 

G(s)  =  -ft}  =  -g-,    -1  „_, i 2_      (J.!) 

V*  +  V-is   +  •••4is  +  qo 

where  the  p    's  and  q   *  s  are  of  the  form: 
*m  ^n 

Pu   =  g^*2   +  ly*    *   iucxy3     *   Ou/3  *  *UJ32  *  lu 

u  =  0,1,2, ...,m  (5-2) 

q     =  a   <*     +  b  cX    +  c  &J6    +  d  8  +  e  Q     +  f 
^•v         v  v  v  /  y  y  v 

v  =  0,1,2, ...,n  (5-3) 


we  may  write: 


m 

£  pusU 


G(s)  =  j^„ 


n 

v 
s 


2   1v 

v=0 


(5-4) 
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Employing  Equation  (5-4)   we  write  in  the  parameterized  form  the 
generalized  compensated  third  order  transfer  function: 


where : 


and: 


^   -    -Iff}"  (5-5) 


P(s)   =  <X2Ks2  +  2<x/3ks  +/S2K  (5-6) 


Q(s)    =  /32[s3  +  £TFr.s2  fZ^r.s  +£TT  r.]  + 

/3[2s4  +  22TTr.s3  +  2STTr.s2  +  2ZTTr.s]  + 

[«5  +  £ipr  .4  +  £TTVS3]  ♦  UTTr.s2  (5-7) 

U  I    1  £    1   -»        3    1 

Making  the  definitions: 

n      i 
A  =  I]   (-l)^bav  ;  etc.  for  Br,Cr,Dr,Er,P  (5-8) 

v=0 
even 

A±  =  £]   (-l)*^"^^  ;  etc.  for  B^C  ,D± ,E.,F.       (5-9) 
v=0 
odd 

m      i 

even 

G.  «  77   (-l)^U-1^J!g  5  etc.  for  H.,I.,J.,K.,L.       (5-11) 
u=0 
odd 

and  substituting  in  Equation  (5-5) >  ^e  write: 

(c*2G     +  K  )      +  j(oy6  I.) 

G(j60h)   =     1  r £ i ^ (5-12) 

(ft\  +/SEr  +  Pp)    +   j(/3\  +/6E.    +  F.) 

Setting  the  magnitude  of  G(j£0,  )  =  M  we  find: 
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,2  .  L,         }I8  .     '^V^'^1/ (5-13) 

b  '        (/3\  +/3Er  ♦  fJ2  ♦  (^\  +/6e±  +  Fl)2 

Manipulating  Equation  (5-13)  algebraically  we  find  that: 

${«,/&)  -£&(*,£)     -  0  (5-14) 
where : 

<j>(<*,/3)  =o<  V  +  2<x2KrGr  +  k2  +  *2/2>2zl  (5-15) 

e(oiyfl)    =/3V   ♦  2/fl3DpEp  *  2yd2DpPp  +/32E2   ♦ 

2/^ErFr  +  *p  +  /^i  +  /^i  +  *i  + 

2i03EiI>i  +  2/92DiFi  +  2^E±Fi  (5-16) 

Substituting  Equations  (5-15)   and  (5-16)    in  Equation  (5-14)   and 
defining: 

P1     =     D2  +  D2  (5-17) 

Qj_     =  2DpEr+  2EJX  (5-18) 

R_      =  2D  F  +  E2  +  E?  +  2D. P.  (5-19) 

R2     =     «2I2  (5-20) 

V,      =  2E  F  +  2B.F.  (5-21) 

1             p.  p     -     i  i  w        / 

¥1     =     ^  +     *i  (5"22) 

W2     =     c<  V     +  2(X2KrGr  +  K2  (5-23) 

we  can  write: 


44. 


m2p.j3?   +  m2qj33    +  fa\  -  R2)/32  +  ifaj3  + 

(M2¥x  -  W2)   =0  (5-24) 

Since  the  Parameter  Pla»e  for  compensation  purposes  has  already 
been  determined  it  is  now  a  matter  of  taking  the  computed  o(.    values 
and  substituting  them  along  with  a  constant  value  of  omega  and  M 
into  Equation  (5-24)  and  then  solving  the  /8  quartic.  This  has  as 
its  solution  the  largest,  real  and  positive  value  of  the  four  roots 
as  before. 

Examples  of  these  curves,  drawn  so  as  to  enable  their  use  as 
overlays  on  Figures  4-2  through  4-10  are  shown  in  Figures  5-1 
through  5-9 • 
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FIG     5-2 


X  -SCALE  =   L00E+01  UNITS'INCH. 
Y -SCALE  =  8B00E+00  UNITS'INCH. 


PRIME,   J„A.,JR.         CONSTANT  BANDWIDTH 
8o0  /  S**3 
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FIG     5-3 

X -SCALE  -   1.S0E+81  UNITS'INCH. 
Y -SCALE  -  2.08E+81  UNITS^INCH. 

PAINE,   JA,JR.         CONSTANT  BANDWIDTH 
125,8  /  S**3       48. 
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FIG     5-4 


X -SCALE  -  9.00E+08  UNITS'INCH. 
Y  -SCALE  -  4.00E+00  UNITS'INCH. 
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FIG     5-5 

X -SCALE  -  4.88E+88  UNITS'INCH. 
Y -SCALE  -  1.00E+88  UNITS'INCH. 

PAINE,  J.A.,JR.    CONSTANT  BANDWIDTH 
■  10/S  ■  »2(S+ 1)     jo. 
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FIG     5-6 


X  -SCALE  =  5.00E+01  UNITS'INCH. 
Y -SCALE  -   1.00E+01  UNITS'INCH. 

PAINE,  JA,  JR.     1/S « «2(S+ 10) 
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FIG      5-7 


K-SCfiLE  -  2.00E+03  UNITS'INCH. 
Y  -SCALE  -  1.00E+02  UMITS/IISCH. 

PAINE,  J0P,?JR0  1/S**2(S+100) 
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FIG  5-8 


X -SCALE  =  3..00E+00  UNITS'INCH. 
Y -SCALE  =  t±n00E+00  UNITS'INCH. 


PAINE,   J.fl.,JR.         CONSTANT  BANDWIDTH 
25,0  /  (S+5.0XS+5.0XS+10.0) 
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fi 


FIG     5-9 


X-SChLE  -  6.09E+99  UNITS'INCH. 


Y -SCALE  -  4.6OE+03  UNITS/INCH. 

PRIME,   J„fl„,JRB         CONSTANT  BANDWIDTH 
25,0  /  (S+5.0XS+5.0XS+ 100.0) 
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6.  Programing  of  the  Compensation  and  Bandwidth  Curves 

Appendices  III  and  IV  are  listings  of  the  programs  developed 
in  the  course  of  investigation  of  this  work.   Only  PROGRAM  PROJECT 
(Appendix  III)  will  be  explained,  as  PROGRAM  PLOT  (Appendix  IV)  is 
essentially  PROGRAM  ROOT  LOCUS  [6]  with  a  minor  subhead  and  internal 
modification. 

PROGRAM  PROJECT  is  designed  to  solve  the  o<  quadratic  and  fi 
quartic  established  in  Sections  4  and  5»  The  program  is  divided 
into  two  main  sections,  the  first  for  the  computation  of  the  o<  -A 
points  and  the  second  for  the  bandwidth  points. 

The  first  section  computes  an  80  by  80  matrix  of  the  Ot  and 
/2    points  corresponding  to  set  values  of  O  and  c*i  ,  The  computa- 
tional part  is  followed  by  two  distinct  graphing  sections,  one  for 
lag  and  the  other  for  lead  compensation. 

The  lag  graphing  section  is  set  up  so  that  during  the  plotting 
of  the  curves  each  value  of  o<  is  tested  to  determine  if  its  value 
is  10  £  °<  ^  1.0001  •   If  no  points  are  found  within  this  range 
then  a  print  out  is  made: 

NO  LAG  COMPENSATION  POSSIBLE 

For  the  lead  section  graphs,  0{   is  again  tested  by  the  criterion 
1.0001^  °^  ^  (X-graph  scale) (X  graph  width)  •  Again  if  there  are  no 
values  of  o(      within  this  region  the  statement: 

HO  LEAL  COMPENSATION  POSSIBLE 
is  printed.   In  this  case  however  a  study  of  the  printed  values  of 

o<   must  be  made  to  insure  that  the  points  are  indeed  non-existant 
or  rather  just  lie  outside  the  range  of  the  graph. 

The  second  main  section  of  the  program  computes  the  value  of 
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/3  for  a  given  value  of  o^  as  determined  by  the  X  graph  scale. 
Here  the  plotting  routine  is  set  up  so  as  to  not  plot  zero  points 
and  to  stop  the  curve  when  either  the  o^    or  /3    value  exceeds  the 
range  of  the  graph. 

There  are  several  go-no  go  options  available  within  the 
program  which  are  amply  explained  in  the  comment  section  preceeding 
the  program  listing  in  Appendix  III, 
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7.  Extensions  to  Higher  Order  Systems 

Although  the  work  presented  to  this  point  has  been  limited  to 
third  order  systems  and  the  program  written  for  this  specific  case, 
investigation  shows  that  generalized  equations  may  be  written  which 
will  allow  the  extension  of  the  program  to  higher  ordered  systems. 
It  can  be  shown  that  for  a  given  n  ~  order  system  with  no  zeros  to 
be  compensated  with  two  identical  sections  of  cascade  compensation, 
that  the  characteristic  equation  of  the  system  may  be  generally 
written  as: 

n+2   „  n+1    2  n   /  2  2   0      2W 

s    +  2ps    +ps  +  (z  s  +  2pzs  +  p  JK  + 

0=n  j=n 

2P  E  (ElTr^  +  p2^   (STTr^  ♦ 

k=n     3  k=n-l     # 

d=i  3=1 

S  (2TTr.)sk  =  0  (7-1) 

k=n+l    £ 
j=l 

where  for  n=4  the  equation  would  be  written: 

6   0  5    24   /  2  2   «       2W 

s  +  2ps^  +  p  s1  +  (z  s  +  2pzs  +  p  )K  + 

2p(E,1!V4  +  2TTVS3  ♦Z'FV2  *ZF")  + 
p'cZTIV3  +£7Tr.s2  +2TTv  +ZTTV)  + 

(STTr.sS  +£^riS4  ♦  Z'JV  +2,5'r.s2)  =  0     (7-2) 

It  may  be  further  shown  that  the  parameters  defined  by  Equations 
(4-4)  through  (4-13)  may  be  written: 
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B 


21 


B 


22 


D, 


D, 


=  KartL(3  ) 

*  Kco2U2("S) 

-2ICcOUo(?) 

-21^1^(5) 


(7-3) 
(7-4) 
(7-5) 
(7-6) 


k=l 


E 


11 


=  2(-i)n+1^n+1un(S)  +  2  f:  [(-i)ko>\^(S )(£Trr.)] 


k=n 
0=1 


(7-7) 


0=n 
k=l 


E 


12 


=  2(.l)n+1^n+\+1(3)   +  2   E[(-l)ka,\(S)(ETrr.)] 


k=n 
0=1 


(7-8) 


0=n 
k=0 


F, 


21 


(-dw.jS)  +  z;  r(-i)kco\1(3)(Zfrri)]+iaj  (3) 

n"x  k=n-l   L  K  x  J        J 

(7-9) 


0=1 


0=n 

k=0 


F, 


22 


=  (-Dna>\(i)  +  e  R-D^VcJxEfrj]  *  kuo(5) 

k=n-l u  J        J 


0=1 


(7-10) 


0=n 


.  (-i)n+2co-\+1( j)  ♦  jg  faf^axzpd 


Go   ■   ( 


k=n+l 
0=1 


k=2 


(7-11) 


-Dn+2o3n+2un+2(S)  ♦  c  [(-D^XcSxEf^j] 


k=n+l 
0=1 


(7-12) 
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These  then  are  the  recursive  equations  required  for  the  complete 

th 
generalization  to  a  n   order  system.  By  employing  the  above 

equations  and  replacing  in  PROGRAM  PROJECT  (Appendix  III)  cards  100 

through  150  and  300  through  540  with  the  appropriate  programing,  the 

program  may  be  used  for  any  given  n   order  system. 

In  like  manner  by  generally  defining: 

P(s)  =  cx2Ks2  +<x/3ks  +/32K  (7-13) 

and:  j=n 

Q(s)    =  sn+2  *  2/3sn+1  +/e2sn  +  2   C    (£TTr.)sk  + 
1  '  k=n         J      1 


d-i 


J=n  o=n 

k=0         ,    k=2 


a2  £  (ETr,)**  +  n  (2TIV)sk  (7-14) 

~  k=n-l   i  k=n+l    3 

and  using  Equations  (5-8)  through  (5-11)  we  may  replace  in  the 

program  cards  2860  and  2880  through  2920,  thus  adapting  this  part 

th 
of  the  program  to  a  general  n   order  application. 
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8.  Conclusions 

Throughout  this  development  of  the  Parameter  Plane  quadratic 
extension,  the  c.  *s  in  the  generalized  coefficient  form: 
n 

Z 

k=0 


£  {\"2  +  V  +  V/*  +  V3  +  V3*  *  sk) 


0 


(8-1) 


have  been  identically  zero.  This  at  first  appearance  might  seem  to 
detract  from  the  generalization.  The  inclusion  of  this  parameter 
does  not  however  introduce  any  great  difficulty  in  the  solution. 
The  change  in  the  development  would  be  to  the  value  of  P,  and  P„ 


which  would  become: 


ci +  fi  h 


P\ 


P2  "  °2  + 
and  the  final  solution  for  c*  which  would  change  to: 

P. 


(8-2) 
(8-3) 


c* 


2B 


21 


pi  -  2B2 A 


4B 


21 


(8-4) 


For  this  case,  new  selection  rules  for  acceptable  values  of  ©< 
would  be  used,  and  would  be  much  like  those  presented  for  Jo   . 

Though  the  extension  of  the  Parameter  Plane  to  include  the 
OL   -  &      quadratic  case  makes  this  tool  even  more  useful,  further 
work  is  still  to  be  done  in  this  field.  Not  only  must  the  equations 
for  the  solutions  of  the  Parameter  Plane  curves  for  such  cases  as: 

»t  ■  V2/82  +  v*/$  +  V*/32  +  v*2  +  V3*  + 

e^/3  *  V<   +  *i/3  *  \  [5]  (8-5) 
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and  higher  ordered  combinations  of  the  parameters  be  developed,  but 
more  efficient  programing  techniques  must  be  developed.   In  the  use 
of  PROGRAM  PROJECT,  for  instance,  as  the  location  of  the  system  poles 
on  the   OT  axis  of  the  S-plane  move  to  the  left,  the  computational 
time  becomes  excessive  due  to  present  programing  technique  and  computer 
speed. 

Another  major  problem  in  further  extensions  of  these  techniques, 
and  indeed  even  other  applications  of  the  curves  from  the  preceeding 
development,  will  be  interpretation.  In  this  case,  the  initial  sub- 
stitution of  variables  immediately  allowed  interpretation  of  the 
curves  sight  unseen.  Here  then,  will  be  most  likely  the  one  single 
drawback  to  further  extension,  for  as  the  parameters  <=*  and  /3     are 
used  as  representations  of  other  variables  in  control  systems,  each 
application  will  have  its  own  unique  interpretation. 
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APPENDIX  I 
THE  ROSS-WARREN  METHOD  OF  COMPENSATION 
The  Ross-Warren  £2j  Method  is  an  algebraic  way  of  determining 
the  location  of  compensator  poles  and  zeros  for  single  or  multiple 
sections  of  cascade  lead  or  lag  compensation.  The  relationships 
of  the  method  are  based  on  relationships  between  the  root  locus  gain 
and  the  error  coefficients,  where  the  uncompensated  root  locus  gain 

±S:  |Hn  TTlr  +  Pil 

|K(r)|   =   ^ (1-1) 

TT  I r  +  z  | 

i=l  ' 

and  the  compensated  root  locus  gain  at  the  point  r  is: 

Kc(r)  =  K(r)  ("IfM  '    («) 

where  k  is  the  number  of  identical  sections  introduced  by  the 

compensator. 

The  other  part  of  the  equality  for  the  derivation  is  determined 

by  the  fact  that  the  uncompensated  error  coefficient  is  defined  by: 

m 

Kx(r)  =  K(r)  -^ (1-3) 

T    p± 

i=l 
and  the  compensated  error  coefficient  is: 

K  (p)-K(r)*— —  (1-4) 

cxv  '  xv  '   z,    p  v   ' 

d    *c 

These  relationships  are  expressed  geometrically  on  the  S-plane  in 

Figure  1-1  and  are  used  to  derive  the  relationships: 


63. 


S-  PLANE       RELATIONSHIPS     TO 
DERIVE      THE      TH£      ROSS-  WARREN    EQUATIONS 


r=-cr,  +jco( 


CO 
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where : 


and: 


Here: 


hK  sin  A  ,    * 

Zc  '   K(r)  sin  0  C3^ 

h   sin  A  ,_  ,s 

pc  °  sin  (8-4.)  (I"6) 


?\  =  cot"   (cot  c£  -  j^j     esc  (()  )  (1-7) 

9  =  Tf  +  <j>  -  J  -  }  (1-8) 


h   -   the  distance  from  the  origin  to  the  root  point 

and  is  specified  by  the  choice  of  zeta  and  omega 

K   -   root  locus  gain  of  the  uncompensated  system 

K(r)  -  root  locus  gain  of  the  uncompensated  system  at 
the  point  r 

CD  -   angle  required  to  make  r  a  point  on  the  root  locus 
/\    -       an  angle  calculable  from  K,  K(r)  and  <j? 
O  -   angle  between  the  negative  real  axis  and  h 
©  -   an  angle  calculable  from  (J>  ,  j  and  /\ 
These  equations  are  generalized  for  n  sections  so  that: 

N  /     lT*        sin  A 

z     =  hY^A-    .         "  (1-9) 

c  1  K(r)  sin  ©  \     * i 

h       sin  A 

*c  ■    Sin(en-^>n)  (I"10) 

where:  n____ 

>n       =      OOt"1    (cot   4>  n   -     Y^7         OBO^n)  (1-11) 

en  -  if  +  4>n-  S  -  >n  d-12) 

and:  , 

dj)  n     =       C|>  /  n  (1-13) 
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APPEiJDIX  II 
DETEBMIMTION  OF  BETA  QUARTIC  SOLUTION 
Here,  three  examples  of  the  approximately  fifty  calculations 
that  were  made  in  the  investigation  of  the  proper  solutions  of  o< 
and  &   ,  will  be  presented.  These  calculations  are  started  by  em- 
ploying Equations  (l-9)  through  (l-13)  with  n  equal  to  two. 
Case  1.     3  =  90°  K  =  30.0    h  =  4.0 

<+>  .  tan-l  1*0       tan"l  ££  +  tan"l  |4  .  180°  ,  24°38» 
T         0.0         1.0         5*0 

<£2  =  12°19« 

K(r)  =  (4.0)(17.0)2(41.0)2  =  105. 60 304 

"Y  */&(*)'   =  .532994 

}\2     =  cot_1(cot  12°19'  -  (. 532994) esc  12019l)  =  25°40t 
62  =  180°  +  12°19'  -  90°  -  25°40»  =  76° 39* 

pc  =^0_Bin2^°40'    .  1#Q225 

sin  64  20' 

2     =    (4.o) (.532994)    &l£5l3Q1       m  #94Q06 

c  sin  76  39' 

this  implies  that: 

/3   =  1.9225  o<   =  2.02568 

Case  2.  "5    =  60°  K-  30.0  h  -  5.0 

<{>    .  ^(J^S)   +  tan     <-4^gt)    +  tan     (-^M)   -  180° 

<fy  =     109°16» 
4>2     =    54°33' 

1  JL 

K(r)   =  (5.0)(21.000199)2(25.00019)a  =  114. 8764 

^(k/kCp)^      =      .511029 
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7\2  =  cot"1  (cot  54°33'  -  (.511029) esc  54°33')  =  85°09.5' 
G2  =  186°  +  54°33'  -  60°  -  85°09.5f  =  89°32.5' 

__    (^o)sm8^.?.   __    8#?2057 

sin  34 50.5 « 


=  (5.0) (.511029)  sin  §S!22±^   =  2#54616 


sin  89°23.5! 
this  implies  that! 

/3  =  8.72057  c*  =  3.42498 

Case  3.       "3  =  90°    K  =  30.0       h  =  5.0 

$  =  *»"*(£§)  +  tan"1^)  +  tan'1^)  -  180°  =  33°48' 

<p2     =  16°54! 

K(r)  =  (5.0)(26.0)2(50.0)2  =  180.277586 
^Iv/K(r)'   =    .407934 

>\2  =  cot"1  (cot  16°54»  -  (.407934) esc  l6°54f)  =  27°59' 
92  =  180°  +  16°54'  -  90°  -  27°59'  =  78°55« 

=  (g-0)   sin  27°^'      =  2#60005 
sin  62  01 » 

an  -  (5.0) (.4079 34)  sin  27n^'    =  .95444 
c  sin  78  55' 

this  implies  that: 

y6  =    2.60005  <*  =    2.72416 

NOW  taking  the  values  of  >0   just  found  we  compare  them  to  the 
roots  of  the  /q     quartic  given  in  Figures  II— lb,  II-2b  and  II— 3b 
respectively.  We  see  that  the  proper  answer  for  /3  >  in  all  cases, 
is  indeed  the  largest  positive  and  real  root  of  the  quartic.   From 


67. 


Figures  II-4  and  II-5  we  now  obtain  the  values  necessary  to  solve 
for  the  value  of  ^  ,  These  steps  for  the  throe  cases  are  shown  in 
Figure  II-6,  and  again  we  note  that  only  when  applying  the  above 
selected  values  of  a3  do  we  obtain  the  correct  value  of  0<  • 
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